Diffusion and transport of spin pulses in an n-type semiconductor quantum well 
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We perform a theoretical investigation on the time evolution of spin pulses in an n-type GaAs 
(001) quantum well with and without external electric field at high temperatures by constructing 
and numerically solving the kinetic spin Bloch equations and the Poisson equation, with the electron- 
phonon, electron-impurity and electron-electron Coulomb scattering explicitly included. The effect 
of the Coulomb scattering, especially the effect of the Coulomb drag on the spin diffusion/transport is 
investigated and it is shown that the spin oscillations and spin polarization reverse along the direction 
of spin diffusion in the absence of the applied magnetic field, which were originally predicted in the 
absence of the Coulomb scattering by Weng and Wu [J. Appl. Phys. 93, 410 (2003)], can sustain 
the Coulomb scattering at high temperatures (~ 200 K). The results obtained are consistent with a 
recent experiment in bulk GaAs but at a very low temperature (4 K) by Crooker and Smith [Phys. 
Rev. Lett. 94, 236601 (2005)]. 

PACS numbers: 72.25.Dc, 72.25.Rb, 73.21. Fg 



Semiconductor spintronics, which aims at replacing the 
charge degrees of freedom with the spin ones in electronic 
devices, has attracted substantial attention recently. 1,2,3 
A thorough understanding of transport of spin polarized 
electrons from one place to another by means of electrical 
or diffusive current is an important prerequisite for the 
realization of spintronic devices such as spin transistors 
and spin valves. Study of the time evolution of a spin 
pulse provides an ideal platform for this purpose and has 
been investigated by two of us based on kinetic spin Bloch 
equations in n-type GaAs quantum wells (QW's). 4,5 It is 
pointed out there that for Zinc-blende semiconductors, 
the spin coherence pker-er = ( c ker c k-<x) pl a Y an impor- 
tant role in the spin diffusion/transport. It leads to some 
striking behaviors, such as at high temperatures (200 K) 
the spin polarization can be opposite to the initial one 
along the spin diffusion even in the absence of an external 
magnetic field. 4 ' 5 This spin oscillation along the direction 
of spin diffusion can not be obtained under the frame- 
work of quasi-independent electron model. 6 This behav- 
ior has been reproduced later by Monte Carlo simulations 
by Pershin. 7 Very recently Crooker and Smith have also 
shown experimentally 8 this feature at very low temper- 
ature (4 K) in an n-type bulk GaAs with strains which 
provide a similar effect as the Dresselhaus term 9 in GaAs 
QW's in our investigation. 

However, it is noted that in our previous spin trans- 
port studies 4,5 as well as the Monte Carlo simulations, 7 
the electron-electron Coulomb scattering is not included. 
Studies of spin kinetics have shown that the Coulomb 
scattering plays a crucial role in spin dephasing and 



relaxation. 10,11,12 ' 13 Therefore, it should also strongly af- 
fect the spin coherence in spin diffusion and transport. 
In addition, the diffusion of a pure spin pulse is formed 
as the spin-up and -down electrons moving in opposite 
directions. This causes inevitably a Coulomb drag 14 ' 15 
between electrons of opposite spins and therefore may 
strongly alter the spin diffusion/transport. All these con- 
siderations inspire us to investigate how the Coulomb 
scattering can affect the spin diffusion/transport. More- 
over, inclusion of the Coulomb scattering further allows 
us to investigate the transport properties of high spin po- 
larization. To our knowledge, a theoretical investigation 
of spin diffusion/transport with the Coulomb scattering 
explicitly included has never been performed in the liter- 
ature. 

We start our investigation in an n-type (001) GaAs 
QW along the z-axis of small well width a. The dom- 
inant spin dephasing mechanism here is the D'yakonov 
and Perel' (DP) mechanism. 16,17 By taking account of 
the DP term, the spin kinetic Bloch equations can be 
written as 4 



<9p(R,k,i) 1 
dt 2 
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with p(R, k, i) standing for the single-particle density 
matrix. The diagonal elements describe the electron dis- 
tribution functions /!w(R, k, t) = /o-(R, k, t) of wave vec- 
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tor k = (k x , k y ) and spin a (= ±1/2) at position R and 
time t. The off-diagonal elements p a _ a (R, k, t) describe 
the inter-spin-band correlations (coherence) for the spin 
coherence. e(R, k,i) = e k S aa ' +h(k) ■cr crcr >/2 — eip(R, t) + 
S CTCT /(R, k, t). Here — k 2 /2m* is the energy spec- 
trum with m* denoting the electron effective mass, er 
are the Pauli matrices. h(k) denotes the effective mag- 
netic field from the DP term which contains contributions 
from both the Dresselhaus term 9 and the Rashba term. 18 
For GaAs QW, the leading term is the Dresselhaus term, 
which can be written as 19 h x (k) = r yk x (k 2 — 7r 2 /a 2 ) and 
h y (k) — ~/k y (ir 2 /a 2 — k 2 ) with 7 standing for the spin- 
orbit coupling strength. 20 The electric potential ip(R,t) 
satisfies the Poisson equation 21 

Vr^(R, t) = e[n(R, t) - no(R)]/(ae ) , (2) 

in which n(R, t) = X] CT k /o-(R-i k, t) is the electron 
density at position R and time t and n (R) rep- 
resents the positive background electric charge den- 
sity. E CTCT < (R, k, t) = - J2q v iPaa> (R, k - q, t) is the 
Coulomb Hartree-Fock self-energy. dp(TL, k, t)/dt\ c and 
dp(R, k, t)/dt\ s in Eq. (1) are the coherent and scattering 
terms respectively. Their expressions are given in detail 
in Ref. 13. For the scattering, we include the electron- 
electron Coulomb scattering, the electron-phonon scat- 
tering and the electron-non-magnetic-impurity scatter- 
ing. For electron-phonon scattering, only the electron- 
longitudinal optical (LO) phonon scattering is consid- 
ered as we concentrate on the high temperature regime 
(200 K) where the electron-acoustic phonon scattering is 
negligible. As we focus on narrow QW's, the separation 
of the subband is large enough so that only the lowest 
subband is needed in our calculation. 

We assume that at initial time t = there is a pure 
spin pulse polarized along z-axis and centered at x = 
0. The electrons are locally in Fermi distribution, i.e., 
f a (x,k,0) = (exp{[e k - fi a (x)]/T e } + 1) _1 , where p, a (x) 
stands for the chemical potential of electrons with spin 
<7 at position x and is determined by the corresponding 
electron density: N rT (x,0) — ^ k f a (x, k, 0). The shape 
of the initial spin pulse is assumed to be Gaussian like 

A7V ff (x,0) = JV V 2 - JV_i/2 - AN e- x2 / Sx2 (3) 

with the A7V and the Sx representing the peak and 
the width of the spin pulse respectively. We further as- 
sume that there is no spin coherence at the initial time, 
p(x,k,0) = 0. 

We numerically solve the kinetic Bloch equations (1) 
together with the Poisson equation (2) following the 
method developed by two of us in Ref. 4. The tricks of 
how to deal with the scattering, especially the Coulomb 
scattering, are laid out in detail in Ref. 13. It is due 
to these tricks which allow us to include the Coulomb 
scattering in the spin diffusion/transport computation. 
By solving the equations, the temporal evolution of the 
electron distribution functions f a (R, k, t) and the spin 
coherence p(R, k, t) at different positions R are obtained 



sclf-consistently. From these quantities, all the transport 
information, such as mobility, particle and spin diffusion 
lengths as well as spin dephasing and relaxation are ob- 
tained explicitly without any fitting parameters and ap- 
proximations such as the relaxation time approximation 
wildly used in the literature. Moreover, as we include 
the Coulomb scattering explicitly, we are able to calcu- 
late the situation with high spin polarization and the sit- 
uation with high external electric field in spin transport, 
where the Coulomb scattering is crucial in the relaxation 
and thermalization of electrons. It is noted that both 
cases have been studied in a much simplified case, i.e., 
pure spin precessions without diffusion/transport (in the 
absence of the spacial degree of freedom). 13,22 

We study the temporal evolution of the spin pulse at 
T = 200 K. The well width and the total electron density 
are chosen to be 7.5 nm and 4 x 10 11 cm -2 respectively 
throughout our calculation. The material parameters of 
GaAs are listed in Ref. 13. The main results of our cal- 
culation are summarized in Figs. 1-4. 




FIG. 1: The absolute value of the spin imbalance |AiV CT | as 
well as the incoherently summed spin coherence p v.s. the 
position x and the time t for the spin pulse with AiVo = 
10 11 cm -2 and 8x = 0.1 /im. Both the electron-LO phonon 
scattering and the electron-electron Coulomb scattering are 
included. 

We first examine if the Coulomb scattering can change 
the spin oscillations along the direction of the spin dif- 
fusion. In Fig. 1 the absolute value of the spin imbal- 
ance \AN a (x,t)\ for the small pulse width of 5x = 0.1 
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/Ltm as well as the incoherently-summed spin coherence 
p{x,t) — ^ k |pkcr-cr(a ; , t)\ are plotted as functions of 
position x along the diffusion direction and the time t 
in the absence of the electric field. Both the electron- 
LO phonon scattering and the electron-electron Coulomb 
scattering are included in the calculation. From the fig- 
ure, it can be seen that the spin signal around the center 
x < 0.06 /im decays very fast due to the diffusion as well 
as due to the spin dephasing. For the region of 0.06 pra 
< x < 0.9 [im, due to the strong diffusion from the center, 
the spin signal first amplifies then decays resulting from 
the weakening of the diffusion as well as the dephasing. 
For the regions further away from the pulse center, i.e., 
0.9 /im< x < 1.02 /im and 2.9 /im< x < 2.94 /mi, we find 
that the above counter effects result in the oscillations of 
the spin polarization with time, which were first discov- 
ered in the absence of the electron-electron scattering. 4,5 
Furthermore it is noted that the striking feature that the 
spin polarization can be opposite to the initial one in the 
absence of a magnetic field is still retained in the presence 
of the Coulomb scattering. It is seen from the figure that 
the spin polarization first reverses at the time about 3 ps 
and the position around 1.02 /zm. And for the region 1.02 
fim < x < 2.96 /xm, one observes a second peak with its 
polarization opposite to the initial one. And then after 
16 ps, the spin polarization reverses back to the initial 
one in the region x > 2.96 [atxi. 

We next investigate the effect of pulse width on the 
spin diffusion by comparing two pure spin pulses with dif- 
ferent pulse widths (Sx = 0.1 /im and 1 /im) but the same 
total spin imbalance (f AN a (x,t = 0)dx = 1.8 x 10 10 
cm -1 ). The spin polarizations at the center of the spin 
pulse is 25% and 2.5% respectively. The spin diffusions of 
the two spin pulses are plotted against t and x in counter 
plots in Fig 2. In the calculation we include the electron- 
electron and electron-LO phonon scattering. It is seen 
from the figure that the decay at the center of the narrow 
spin pulse is much faster than that of the wide one and 
that it is much faster for the narrow spin signal to diffuse 
out. Also the diffusion of the narrow spin pulse shows 
stronger spin oscillations along the direction of the diffu- 
sion. The reverse of the spin polarization occurs far away 
from the position of the initial spin pulse. Moreover, the 
spin polarization reverses back to the initial direction for 
the narrow spin pulse at positions farther away from the 
initial pulse. This peak can hardly be seen for the case of 
the wide pulse. These features are understood as follow- 
ing: The narrower spin pulse leads to larger gradients of 
the spin densities and hence a larger broadening in the 
/c-space, and then leads to a faster diffusion. Moreover, 
a larger momentum k indicates a stronger effective mag- 
netic field from the DP term, which results in a stronger 
spin oscillations. 

In Fig. 3 we investigate the effect of the Coulomb 
drag 14,15 on the spin diffusion by plotting the absolute 
values of the spin imbalance as functions of x at t = 5, 15, 
35 and 60 ps, with and without the Coulomb scattering. 
Contrary to the common impression that the Coulomb 
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FIG. 2: (Color online) The contour plots of spin imbalance 
AN v.s. the position x and the time t for different pulse 
widths, (a): Sx = 0.1 fim; (6): Sx = 1.0 fim. T = 200 K. 



scattering conserves the total momentum and should not 
change the spin diffusion dramatically, one finds that 
when the Coulomb scattering is included, both the dif- 
fusion of the spin pulse and the decay of the spin po- 
larization at the center of the pulse become much slower. 
Specifically at t = 15 ps, a third peak with the spin polar- 
ization along the initial one has appeared when only the 
electron-LO phonon scattering is included. This peak 
is absent when the Coulomb scattering is added. At 
t = 35 ps, the third peak also appears in the case with the 
Coulomb scattering, but with a smaller magnitude. Later 
at t = 60 ps, it is interesting to see that the magnitude 
of the third peak in the case with the coulomb scattering 
exceeds the one without, which indicates a longer spin 
diffusion length with the Coulomb scattering. 

These features are understood as following: The inter- 
change of the momentums between the oppositely moving 
spin-up and -down electrons due to the Coulomb scatter- 
ing induces an effective friction, which is known as the 
Coulomb drag. 14,15 This drag weakens the diffusion of 
the spin pulse dramatically as shown in Fig. 3. In addi- 
tion to the Coulomb drag, the Coulomb scattering also 
markedly affects the spin dephasing. 12 It drives electrons 
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3: (Color online) The absolute values of the spin imbal- 
ance |AiVo-| v.s. the position x at certain times t for the cases 
with only the electron-LO phonon scattering (orange); the 
electron-phonon and electron-electron scattering (blue); and 
the electron-phonon, electron-electron and electron-impurity 
scattering (green), the red curve shows the initial spin pulse. 
AN = 10 11 cm~ 2 , 8x = 0.1 ^m and Ni = 0.5N e . Solid 
curves: AiV > 0; Dashed curves: A7V < 0. 



to a more homogeneous states in the k-space, suppresses 
the anisotropy caused by the DP term, and reduces the 
spin dephasing. 12 Both effects slow down the decay of the 
spin signal at the center of the pulse. Moreover, the spin 
diffusion length is determined by the above two counter 
effects: On one hand, the Coulomb drag shortens the spin 
diffusion length; On the other hand, the increase of the 
spin dephasing time by the Coulomb scattering leads to 
a longer spin diffusion length. And our results indicate 
that for the case of our study, the later effect is more 
important in the long run, which gives rise to a longer 
spin diffusion length. 

We further explore the effect of impurities on spin dif- 
fusion in the presence of the Coulomb scattering. The 
result is also plotted in Fig. 3. From the figure one can 
see that as expected, the spin diffusion becomes much 
slower due to the decrease of the mobility. This effect, 
together with the well known fact that the spin dephasing 
time is increased by the electron-impurity scattering, 17 ' 23 
make the decay of the spin polarization at the center of 
the pulse much slower. 

Finally we investigate the transport of the spin pulse 
under an external electric field E = 500 V/cm along the 
— x-direction. The hot-electron effect is not important 
under this field. 13 Figure 4 shows the differences of the 
spin-up and -down electron densities versus the diffusion 
length x at t = 6, 12, 16 and 24 ps in different cases such 
as with/without the Coulomb and impurity scattering. 
From the figure, one finds that due to the electric field, 



the spin signals do not diffuse symmetrically around the 
center any more, but transport against the direction of 
the electric field. Moreover, the features of the trans- 
port /diffusion of the spin polarization in the side against 
the external electric field become much richer than the 
side parallel to the field. All these results are consistent 
with what observed in the experiment. 8 From the figure, 
one can further see the effect of the Coulomb scattering 
by comparing the curves with (solid curve) and without 
(dashed curve) the Coulomb scattering. It is found that 
similar to the field free case, the spin signal is more dif- 
ficult to diffuse/ transport out due to the Coulomb drag. 
Moreover as the Coulomb scattering impedes the spin 
dephasing, combined with the drag effect, it leads to the 
slower decay of the spin signal at the center. When the 
impurity with impurity density Ni = 0.5iV e is added to 
the system, it further slows the spin transport and the 
decay of the center of the spin signal as the impurity scat- 
tering impedes both the mobility and the spin dephasing. 
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FIG. 4: The spin imbalance AN a v.s. the position x at certain 
times t. Dotted curve: the initial spin pulse; Dashed curve: 
with only electron-LO phonon scattering; Solid curve: with 
electron-LO phonon and electron-electron Coulomb scatter- 
ing; Chain curve: with electron-LO phonon, electron-impurity 
and electron-electron Coulomb scattering. AiVo = 10 11 cm -2 , 
Sx = 0.1 jLtm and Ni = 0.5iV e . 

It is noted that in our two-dimensional case, the tem- 
perature is much higher than the one in the experiment 
(4 K). 8 It is because the system in the experiment is a 
bulk one where the effective magnetic field from the DP 
term and the strain is as weak as ~ 1.6 x 10~ 3 T. It 
therefore requires a much lower temperature to reduce 
the scattering in order to observe the reverse of the spin 
polarizations. While in QW's with small well width, due 
to the strong confinement, the effective magnetic field is 
about 3.4 T and the phenomenon can sustain a much 
stronger scattering at high temperatures. 
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In conclusion, we have performed a study on the 
time evolution of spin pulses at a high temperature 
(200 K) through self-consistently solving the kinetic spin 
Bloch equations with the electron-LO phonon, electron- 
impurity and electron-electron Coulomb scattering ex- 
plicitly included. The effect of the Coulomb scattering 
on the spin diffusion/ transport is clearly demonstrated: 
The Coulomb scattering on one hand prolongs the spin 
dephasing time, which helps to sustain the spin signal. 
On the other hand, it causes the Coulomb drag which 
weakens the diffusion of the spin signal. In the case of 
our study, the Coulomb scattering leads to a slower spin 
diffusion speed but a longer spin diffusion length. 

It is also shown that the Coulomb scattering does not 
kill the striking feature of the spin oscillation and the 
spin reverse along the spin diffusion direction, which were 
originally predicted where only the electron-LO phonon 
scattering is considered. 4 In addition, when an external 
electric field is applied, the spin signal tends to trans- 
port against the direction of the electric field. All these 



phenomena have been observed in a recent experiment in 
bulk GaAs at very low temperature. 8 Our study shows 
that one may observe these features at much higher tem- 
perature in QW's with small well width. Experiments 
are suggested to verify these effects. 

Moreover we also study the effects of the pulse width 
and the impurity on the spin diffusion. The results in- 
dicate that the wider the pulse width is, the weaker the 
spin diffusion and the fewer the times of the reverse of 
the spin polarization along the diffusion. Also the impu- 
rity affects the spin diffusion in two competing directions, 
i.e. impeding the spin diffusion and reducing the spin de- 
phasing. 
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